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Electrical contacts to semiconductors play a key role in electronics. For nanoscale electronic 
devices, particularly those employing novel low-dimensionality materials, contacts are expected to 
play an even more important role. Here we show that for quasi-one-dimensional structures such 
as nanotubes and nanowires, side contact with the metal only leads to weak band re-alignement, 
in contrast to bulk metal-semiconductor contacts. Schottky barriers are much reduced compared 
with the bulk limit, and should facilitate the formation of good contacts. However, the conventional 
strategy of heavily doping the semiconductor to obtain ohmic contacts breaks down as the nanowire 
diameter is reduced. The issue of Fermi level pinning is also discussed, and it is demonstrated 
that the unique density of states of quasi-one-dimensional structures make them less sensitive to 
this effect. Our results agree with recent experimental work, and should apply to a broad range of 
quasi-one-dimensional materials. 



The early work of Schottky, Mott and Bardeen has laid 
the course for much of the fundamental understanding 
and improvement in the performance of electrical con- 
tacts to bulk semiconductors. For nanoelectronics, con- 
tacts are a significant fraction of the total device size and 
are expected to be crucial for device behavior. However, 
because of the unique properties of nanostructures, much 
of the basic fundamental aspects of contacts need to be 
re-examined at the nanoscale. 

The current understanding of contacts to nanostruc- 
tures is in its infancy, both from an experimental and 
theoretical perspective. An example is carbon nanotubes 
(NTs): despite much experimental work, it is still unclear 
wether the contacts are Schottky or ohmic, with reports 
of Schottky contacts for Ti [1] and ohmic contacts for Au 
[2] and Pd [3,4]. However, recent experimental work [5,6] 
has suggested that the type of contact depends on the 
NT, with Schottky contacts for small diameter NTs and 
ohmic contacs for large diameter NTs. 

From a theoretical perspective, it has been demon- 
strated that the concept of Fermi level pinning (crucial 
in traditional semiconductors) is ineffective for quasi-one- 
dimensional nanostructures end-bonded to metals [7]. For 
NTs side- contacted by a metal, modeling has been used 
to extract Schottky barriers from experimental measure- 
ments of the ON current in NT transistors [5] , but have 
not addressed the origin of the Schottky barriers; and 
atomistic calculations have provided case-by-case stud- 
ies [8,9] for planar contacts. However, a more general 
theoretical understanding for side-contacts to quasi-one- 
dimensional (QID) structures is still missing, especially 
in light of the recent experimental findings. 

In this paper, we present a theoretical and modeling 
analysis of side contacts to nanotubes and nanowires. We 
show that the conventional concepts developed for bulk 
metal-semiconductor contacts do not simply carry over 
to the nanoscale. In particular, band re-alignement due 
to charge transfer is weak due to the limited available 
depletion width. In NTs, this leads to relatively small 



and slowly varying Schottly barriers with NT diameter. 
In nanowires (NWs), there is a range of diameters with 
minimized Schottky barriers, providing optimal contact 
properties. We also demonstrate that in general, QID 
structures are much less sensitive to Fermi level pinning 
than their bulk counterparts. Finally, a conventional 
strategy for making ohmic contacts is to heavily dope the 
semiconductor near the contact; we show that at typical 
dopings for Si, the contact resistance increases rapidly as 
the nanowire diameter is decreased below 40 nm. 
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FIG. 1. Panel (a) shows a cross-section of the 
metal-nanostructure contact along the length of the nanos- 
tructure. Panels (b) and (c) show radial cross-sections for 
metal-nanotube and metal-nanowire contacts. The NT is 
separated from the metal by a distance s = 0.3 nm. The 
NW is modeled as a solid cylinder with a sharp interface 
with the metal. Panel (d) shows the band alignment at a 
metal-semiconductor contact before charge transfer. In a bulk 
contact, panel (e), band bending over a distance W due to 
charge transfer leads to a Schottky barrier At^jfe. 

We begin by describing the contact geometry consid- 
ered here. Figure la shows a sketch of a cross section 
of the contact consisting of a QID structure embedded 
in a metal. For explicit systems, we consider a single- 
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wall NT , as shown in Fig. lb, or a solid nanowire as in 
Fig. Ic. For the NT, the metal forms a cylindrical cavity 
of radius R + s where R is the NT radius and s is the 
distance between the NT and the metal, while for the 
NW we consider a solid, continuum cylinder embedded 
in a perfect metal, with a sharp interface between the 
nanowire surface and the metal. 

In the simplest picture, the difference between the 
metal Fermi level Ep and the semiconductor valence 
band edge Ey (the barrier for holes) is simply given by 
(Fig. Id) 



Ao = Ef - Ey = -E„ + 



(1) 



where $m and $s are the metal and semiconductor work- 
functions respectively, and Eg is the semiconductor band 
gap. A positive value for Aq indicates a Schottky bar- 
rier, while a negative value indicates an ohmic contact. 
Because bandgap decreases with increasing diameter for 
QID structures, the value of Aq depends on the nanos- 
tructure diameter. The behavior of Eq. (1) for undoped 
NTs is shown in Fig. 2 for a value of — (f)^y = 0.4 eV 
(typical of Pd), and using the relation Eg = 2a'y/d be- 
tween bandgap and NT diameter d {a = 0.142 nm is the 
C-C bond length and 7 — 2.5 eV is the tight-binding over- 
lap integral). One problem with this picture (besides the 
fact that the physics is incomplete, as will be discussed 
below) is that Eq. (1) predicts large and negative values 
for Ao, signaling strong ohmic contacts. However, it is 
clear that such strong ohmic contacts are not observed 
experimentally. 

In general, because the metal Fermi levels of the metal 
and semiconductor are not equal, charge transfer between 
the metal and semiconductor occurs, and leads to band 
re-alignement. At a bulk semiconductor junction (Fig. 
le) this charge transfer leads to the Schottky barrier 



^bulk 



E„ 



(2) 



where x is the semiconductor electron affinity. This re- 
lationship arises because, in the bulk system, a deple- 
tion width perpendicular to the metal-semiconductor in- 
terface is created in the semiconductor until the band 
lineup in Eq. (2) is obtained. However for QID struc- 
tures, the depletion width depends exponentially on the 
doping [10] and is much longer than the device size for 
non-degenerate doping, leading to slowly varying bands 
outside of the contact; and for a three-terminal device 
the band-bending in the channel is governed by the gate 
voltage. In either case, the band alignment is deter- 
mined by that in the contact. But for a side-contacted 
QID structure, the semiconductor is only a few nanome- 
ters thick in the direction perpendicular to the metal- 
semiconductor interface; thus only a region of the order 
of the nanostructure cross-section can be depleted, giving 
partial band re-alignement. The value of A will then be 



somewhere between Aq and Ahuik (for an undoped NT or 
NW, 

Abuik — which would always give relatively 

high Schottky barriers). 

Nanotubes are an extreme example of this situation, 
since the possible "depletion width" is the size of the NT 
wall; the charge transfer and image charge in the metal 
create two nested hollow cylinders with opposite charge, 
and an associated electrostatic potential. This electro- 
static potential in turn shifts the bands, and changes the 
amount of transferred charge. Thus, the charge and po- 
tential must be determined self-consistently. We can cap- 
ture this behavior using analytical models for the charge 
and potential. The charge per unit area on the NT can 
be expressed as 



a = eN j DNT{E)f{E - EF)dE 



where 



Dnt{E) 



a\/3 



\E + eVNT\ 



^'^yiE + eV^rf-iEg/2f 



(3) 



(4) 



is the NT density of states [11], f {E — Ep) is the Fermi 
function, and N = 4/(3\/3a^) is the atomic areal den- 
sity. We assume a uniform and sharp distribution of the 
charge on the NT. 

For the geometry of Fig. 1, solution of Poisson's equa- 
tion gives the potential on the NT as 



eR^ R + s 
bVnt = -cr — m — — = -770- 

£0 -K O 



(5) 



where £0 is the permittivity of free space and C is the ca- 
pacitance per unit area between the metal and the NT. 
Equations (3) and (5) can be solved self-consistently for 
a given NT. In this model the electrostatic potential in- 
duced on the NT modifies the barrier to 



A = Aq — bVnt- 



(6) 



Figure 2 shows results of such calculations for parameters 
typical of Pd. Clearly, the behavior is different from the 
simple expressions in Eqs (1) and (2). The results suggest 
that there is a transition between Schottky and ohmic 
behavior at a NT diameter around 1.4 nm, in agreement 
with recently published experimental data for Pd con- 
tacts [5,6]. 

The results of these calculations can be vcirified using 
an atomistic treatment of the NT. For selected zigzag 
NTs, we apply a tight-binding Green's function formal- 
ism, dividing the system into semi-infinite left and right 
leads connecting a central scattering region (Fig. la). In 
the central region, we calculate the Green's function from 



= [{E -eV)I-Ho- E«] ' 



(7) 



where Hq is the tight-binding Hamiltonian for the iso- 
lated NT and V is the electrostatic potential on each 
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ring of the zigzag NT. The function represents the 
interaction of the scattering region with the semi-infinite 
NT leads. To obtain the charge density, we assume a 
uniform distribution of the charge in the azimuthal di- 
rection, and spatially distribute the total charge on each 
ring in the radial and axial directions with a Gaussion 
smearing function. Thus the three-dimensional charge 
density is given by 



p(r, 0, z) 



I 



g{z- zi,r- 



R)^j dEImGjj (8) 



where g(x,y) = {An'^RazCrr) ^ exp [-x'^ /2al — y'^ /2a'^] 
with R the tube radius, zi the position of ring I, and cr^ 
and ar the smearing lengths in the axial and radial di- 
rections respectively (this expression for g is valid when 
R^ ar, and we use values of cr^ = 0.14 nm and ct^ = 0.06 
nm). 

The electrostatic potential is calculated by solving 
Poisson's equation in cylindrical coordinates on a grid 
with the geometry of Fig. 1. Once the three-dimensional 
electrostatic potential is obtained, the value for V on 
each layer is taken as the value of the three-dimensional 
electrostatic potential at the atomic position of each ring 
along the NT. We iterate the calculation of the charge 
and the electrostatic potential until self-consistency. 

As shown in Fig. 2, the results of such calculations in- 
dicate excellent agreement with the analytical approach 
introduced above. This agreement allows us to use equa- 
tions (3) and (5) to gain further understanding of the 
properties of these contacts, as we now discuss. 
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2. Schottky barrier at nanotube-metal contacts for 



parameters typical of Pd. Dotted line is prediction from Eq. 
(1), solid line is computed from Eqs (3) and (5), dashed line is 
Eq. (11) and circles are calculated from atomistic approach. 

To proceed further we focus on the small and posi- 
tive A regime, and approximate the fermi distribution as 
J{E~Ef) w \eyip[-l3{E - Ef) /kT] where /3 = 0.7. 
Asymptotic expansion of the integral in Eq. 3 leads to 
the charge 



eiVaVS lEgkT _^ 



e '^''T . 



(9) 



Combined with Eq. (5) this expression for the charge 
gives the Schottky barrier 



(10) 



where a = {e^NaVS) / {2y/2^TT^/^R'yC) and L{x) is 
Lambert's W function. A more appealing formula can 
be obtained using asymptotic expansion of L{x) giving 



(11) 



The behavior of this function is plotted in Fig. 2, showing 
good agreement with the full calculation. The logarith- 
mic dependence implies relatively slowly varying A, at 
least compared with Eq. (1). The NT diameter delimit- 
ing Schottky from ohmic behavior is [12] 




1 + a 



kT 



NT 



(12) 



thus the diameter is increased from its bare value by 
5d = Q^^/^~^~7'^o- Making ohmic contact to a wide 
range of NT diameters requires a small 5d] this can be 
accomplished at low temperature, or with a large metal 
workfunction. Because a is inversely proportional to the 
capacitance, our result also provides one reason why em- 
bedding the NT in the metal and heating can improve 
the contact properties: the embedded contact provides 
a larger capacitance compared to a planar contact (for 
typical parameters we find a « 2.5 for the embedded 
contact and up to a w 14 for a planar contact), while 
heating presumably provides a smaller average value for 
s and a smaller capacitance. The large value of a for 
planar contacts imply that this type of contact is almost 
pathological, since the value of 5d should be very large 
unless the metal worfunction itself is large. This may 
explain why ab initio calculations find Schottky barriers 
[8,9] despite the large values of — ^mt used in the 
calculations. It is worth mentioning that changing the 
contact geometry to modify the contact properties is not 
possible in bulk contacts, and is thus a unique feature of 
nanostructures . 

We now consider side-contacts to nanowires. The dif- 
ference between NWs and NTs is two-fold: first, in NWs 
the possible depletion width increases with diameter, 
while in a NT it stays constant; second, while the band 
gap also decreases with diameter for NWs [13], the large 
diameter limit can lead to a finite band gap, while for 
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NTs it leads to a zero band gap. It is not clear a priori 
how these effects influence the contact behavior. To ad- 
dress this issue, we consider a model NW with density of 
states 



Dnw{E) 



Tlfl 



(E - EJ2) 



-1/2 



(13) 



where m* is the effective mass. For silicon NWs, it has 
been shown experimentally [13] that the band gap de- 
pends on diameter as Eg = Eq + C / d? where Eq 1.12 
eV and C = 4.33 eVnm^. We consider the situation 
— ^NW — 0.7 eV typical of contacts to Si. Fig. 3a 
shows the expected Schottky barrier heights from Eq. 
(1) , which predicts ohmic contacts to NWs with diam- 
eters larger than 4 nm. To study the effects of charge 
transfer, we perform a self-consistent calculation of the 
charge and potential, using Eq. (13) to obtain the charge 
and solving Poisson's equation numerically in the NW to 
obtain the potential (we use an atomic volume density 
iV„ = 5 X 10^* atoms/m^). 
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FIG. 3. Panel (a): Schottky barrier at nanowire-metal con- 
tacts for parameters typical of SiNWs. Dotted line is predic- 
tion from Eq. (1) and solid line is self-consistent calculation. 
Panels (b) and (c): Band-bending across nanowires with di- 
ameters of 2 nm and 10 nm, respectively; dotted lines are 
Fermi level. 

Fig. 3b, c shows the self-consistent band-bending for 
NWs of 2 and 10 nm radius. Clearly, the nanoscale 
dimension of the NWs prevents the bands from reach- 
ing their asymptotic value; instead, there is only a rel- 
atively weak band-bending present. To quantify the 
Schottky barrier height, we calculate the spatial aver- 



age of Ep — Ey{r)] the results plotted in Fig. 3a in- 
dicate a remarkable behavior: the contact is always of 
Schottky character, with the barrier minimized at a di- 
ameter of about 4 nm. Thus, while in NTs the barrier 
height decreases monotonically with diameter, the behav- 
ior in other QID structures may be non-monotonic, with 
a range of diameters providing optimal contact proper- 
ties. 

We now discuss the issue of Fermi level pinning. As 
shown in Fig. 4a, in a bulk metal-semiconductor contact, 
metal-induced gap states (MIGS) lead to band-bending 
over a distance I <C VF, and modification of the Schottky 
barrier height to Api„. The question is to what extent 
this mechanism influences the properties of side contacts 
to NTs and NWs. To model this effect, we consider a 
radial pinning charge [7] 



^ pin 



{r)=DoNA[EF~EN{r)]h{r) 



(14) 



where the neutrality level En is at midgap [i.e. Eiy{r) — 
-eV{r)], h{r) = e"''/' for a NW and /i(r) = 6rM for 
a NT, and Na = N for a NT and Na = N^^^ for 
a NW. We choose I = 0.3 nm, a typical value for 
metal-semiconductor interfaces [14]. We add this pin- 
ning charge to Eq. (3) or to the charge calculated from 
Eq. (13) and repeat our self-consistent calculations. 

Figure 4b shows the Schottky barrier calculated for 
several NTs as a function of the density of gap states 
{Apin = Eg/2). Clearly, there is a rapid onset of pin- 
ning at Dq ~ 0.1 states/(atom-eV); this value of Dq 
is rather larger considering the van der Waals bonding 
of NTs to metals. Furthermore, atomistic calculations 
[8] have obtained values at least an order of magnitude 
smaller. Thus, as in end-bonded contacts, we expect that 
Fermi level pinning will play a minor role in side-contacts 
to NTs. The rather large value of Dq required to see pin- 
ning effect is a generic property of QID structures as we 
now discuss. 

Figure 4c shows the effects of Fermi level pinning 
on the barrier height in SiNWs. The results also in- 
dicate a value of Dq ^ 0.1 states/(atom-eV) required 
to see pinning effects. For comparison, the inset in 
this flgure shows the same calculation for a bulk metal- 
semiconductor interface with the same parameters, in- 
dicating that only 0.002 states/ (atom-eV) are needed to 
reach the onset of pinning. Thus, the QID system re- 
quires almost two orders of magnitude larger density of 
pinning states compared with the bulk interface. 

The origin of this behavior can be traced to the unique 
density of states of QID systems. Indeed, for Si, we 
can repeat the analysis leading to Eq. (9) using the 
density of states for the NW and for the bulk system 
[DbuikiE) = V2im*f^ i^^^Y' - Eg/2]. This 
leads to the ratio aNw/<^buik — (2t:Nv^^P^ / {m*kT). 
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FIG. 4. Panel (a) shows a sketch of the band-bending 
in the presence of metal-induced gap states at a bulk 
metal-semiconductor interface. Panels (b) and (c) show the 
calculated Schottky barrier as a function of the density of gap 
states for several NTs and NWs, respectively. The inset in (c) 
shows the behavior for a planar metal-semiconductor contact. 



The appearance of the kT factor in the denominator is 
entirely due to the QID density of states of the NW and 
the presence of a van Hove singularity at the band edge. 
At room temperature, we find that (Tjviy/cbuifc > 100; 
thus the MIGS are competing with a much larger charge 
density in the QID system. 
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FIG. 5. Band-bending across Si NWs with doping of 
10^^ e/cm~"^ for diameters of 40 nm (a) and 10 nm (b). The 
arrow indicates tunneling of electrons throught the Schottky 
barrier. The normalized resistance is shown in figure (c) as a 
function of NW diameter. 

Our discussion has so far focused on the situation of 
low doping, where the strategy for making ohmic contacts 
is by selection of a metal with appropriate workfunction. 
In traditional metal-semiconductor contacts, an alterna- 
tive approach is to heavily dope the semiconductor, and 
rely on tunneling through the Schottky barrier to reduce 
the contact resistance and obtain ohmic-like contatcs. To 
address the feasibility of this approach for contacts to 
NWs, we repeat our self-consistent calculations for the 
Si NW, focusing on the situation where the metal Fermi 
level is in the middle of the NW bandgap at the inter- 
face, and adding a uniform doping charge of 10^^e/cm~^. 
Figure 5 shows the band-bending in the presence of this 
doping charge for NWs of 40 and 10 nm diameters. We 
calculate the contact conductance from 



G 



dE 



(15) 



where the tunneling probability T(E) is obtained from 
the WKB approximation. The normalized contact resis- 
tance is then Goo/G where Goo is the conductance in the 
limit of large diameters. The behavior of the normalized 
resistance as a function of NW diameter is shown in Fig. 
5c, indicating a rapid increase of the resistance with de- 
crease in diameter. The origin of this behavior is due 
to the increased tunneling distance and reduced range of 
tunneling energies because of the poor band-bending in 
the NW (the increase in bandgap also plays a role for 
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the smaller diameters). One implication of this result is 
that different diameter NWs will require different doping 
levels to achieve the same contact quality. 

In summary, wc find that the concepts developed to de- 
scribe traditional metal-semiconductor interfaces fail to 
properly account for the properties of contacts to quasi- 
one-dimensional structures such as nanowircs and nan- 
otubes. The nanometer cross-section of nanowires and 
nanotubes prevents the formation of an appropriate de- 
pletion width, leading to only weak band re-alignemcnt; 
and because of the diverging density of states at the 
band edge, these nanostructures are much less sensitive 
to Fermi level pinning than their bulk counterparts. Op- 
timizing device performance will not only require select- 
ing QID structures for their behavior in the channel, but 
also for their contact properties. We expect that our 
results will be applicable to a broad range of QID struc- 
tures. 
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